Copy 


NASA MEMO 10-23 -58L 




NASA 


MEMORANDUM 


EFFECTS OF BODY SHAPE ON THE DRAG OF A 45° 

SWEPT BACK -WING— BODY CONFIGURATION AT 

MACH NUMBERS FROM 0.90 TO 1.43 

By Walter B. Olstad and Thomas L. Fischetti 

Langley Research Center 
Langley Field, Va. 


NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 


WASHINGTON 

November 1958 

Declassified April 12, 1961 




fill"! Ill 






NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


NASA MEMO 10-23-58L 


EFFECTS OF BODY SHAPE ON THE DRAG OF A 45° 
SWEPTBACK-WING — BODY CONFIGURATION AT 
MACH NUMBERS FROM 0-90 TO 1.43* 

By Walter B. Olstad and Thomas L. Fischetti 

SUMMARY 


An investigation was made of the effects of body shape on the drag 
of a 45° sweptback-wing — body combination at Mach numbers from 0.90 

to 1.43. 

Both the expansion and compression fields induced by body indenta- 
tion were swept back as the stream Mach number increased from 0.94. The 
line of zero pressure change was generally tangent to the Mach lines 
associated with the local velocities over the wing and body . The strength 
of the induced pressure fields over the wing were attenuated with span- 
wise distance and the major effects were limited to the inboard 60 per- 
cent of the wing semi span. 

Asymmetrical body indentation tended to increase the lift on the 
forward portion of the wing and reduce the lift on the rearward portion. 
This redistribution of lift had a favorable effect on the wave drag due 
to lift. Symmetrical body indentation reduced the drag loading near the 
wing-body juncture at all Mach numbers. The reduction in drag loading 
increased in spanwise extent as the Mach number increased and the line 
of zero induced pressure became more nearly alined with the line of 
maximum wing thickness. 

Calculations of the wave drag due to thickness, the wave drag due 
to lift, and the vortex drag of the basic and symmetrical M 3 1.2 body 
and wing combinations at an angle of attack of 0° predicted the effects 
of indentation within 11 percent of the wing— basic-body drag through- 
out the Mach number range from 1.0 to 1.43- Calculations of the wave 
drag due to thickness, the wave drag due to lift, and the vortex drag 
for the basic, symmetrical M = 1.2, and asymmetrical M = 1.4 body and 
wing combinations predicted the total pressure drag to within 8 percent 


* 
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of the experimental value at M = 1.43. The incremental pressure drag 
due to a symmetrical M = 1.2 indentation wan predicted to within 8 per- 
cent of the total pressure drag of the wing-indented body combination for 
all Mach numbers above M = 1.0. 

A theoretical method presented herein proved useful in predicting 
the distributions of the oblique section lift parameter ; |3Z/2q. The 

differences in wave drag due to lift calculated from the predicted and 
experimental distributions of 3 i/ 2q, were lees than 6 percent of the 
total pressure drag. 


INTRODUCTION 


The transonic area rule (ref. l) provides a means for reducing the 
zero-lift wave drag of wing-body combinations at Mach numbers near 1. 

An extension of the area-rule concept led to a supersonic area rule 
(refs. 2 and 3) that has proven useful in reducing the zero-lift wave 
drag at higher speeds. The zero- lift wave drag of a wing-body combina- 
tion is related, by the supersonic area rule, to the wave drags of a 
series of equivalent bodies. The equivalent-body shapes are determined 
by the area distribution of the wing-body combination. The area rule 
as presented in references 2 and 3, neglects the reflection of disturb- 
ances by the wing, but for unsymmetrical configurations, this intro- 
duces errors. Reference 4 suggests including these effects by con- 
sidering separately the area distributions above and below the wing- 
chord plane. Application of this concept decreased the drag signifi- 
cantly for two unsymmetrical configurations. 

Reference 5 considers the wave drag of a Lifting wing-body combina- 
tion at supersonic speeds. The wave drag is again related to the wave 
drags of a series of equivalent bodies of revolution, but in this case 
the equivalent-body shapes depend on both the pressure distribution and 
the area distribution of the wing-body combination. 

In the past, experimental studies of the area rule have been con- 
fined largely to measurements of overall force.;. The present paper dis- 
cusses the effects of body shaping on the distribution of drag. A 45° 
sweptback wing has been tested in combination with four different bodies 
at a Mach number of 1.43 and with two of these bodies at Mach numbers 
from 0.8 to 1.43. The pressure distributions ?ire studied to determine 
the magnitude and extent of the pressure field due to both symmetrical 
and asymmetrical body indentations and to determine the locations on 
the wing and body where the drag is reduced as a result of the indenta- 
tions. In addition, experimental drag values for asymmetrical and lifting 
configurations are compared for the first time with values calculated by 
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methods based on references 4 and 5 and by a modified method presented 
herein. 


SYMBOLS 


b 

c 

c 

c av 

c d 

C D 

c D,f 

C D,p 

c D,w 0 

C D,w 

c D,w,6 

C D ,v 

c L 



h 


wing span, 28.478 in. 
local chord, in. 

mean aerodynamic chord, 8,42 in. 
average wing chord, 7-12 in. 

section drag coefficient 

Draff 

total drag coefficient, m 

skin- friction drag coefficient 
integrated pressure drag coefficient 
wave drag coefficient due to thickness 
wave drag coefficient due to lift 
incremental wave drag coefficient due to lift 

vortex drag coefficient 

T i ft 

lift coefficient, Sw - 
section normal-force coefficient 
local pressure coefficient 
differential pressure coefficient 
section chord-force coefficient 
vertical dimension, in. (see fig- 8) 


k 


l or l(x,6) oblique section lift 
M free- stream Mach number 

<1 free- stream dynamic pressure 

R body radius 


R 


max 


maximum body radius 



0 cutting-plane roll angle, deg 

A angle between the wing leading edge ind the y-axis, deg 

1 auxiliary coordinate in x-direction, in. 

t angle between intersection of Mach cutting plane with hori- 

zontal plane and x-axis 

3Z/2q oblique-section lift parameter, in. 

Subscripts : 

a due to asymmetry about wing- chord plane 
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B body 

t due to thickness 

W wing 

WB wing-body combination 

a due to angle of attack 

Primes indicate differentiation with respect to x. 


MODEL CONFIGURATIONS 


The wing used in this investigation had an aspect ratio of 4.0 and 
a taper ratio of 0.15. The wing section varied linearly in thickness 
from a NACA 65A206, a = 0 section at the root to a NACA 65A203, 
a = 0.8 (modified) section at 50 percent of the semispan. Then the 
wing section remained constant to the tip. The wing was tested in con- 
junction with four body shapes. A sketch of the wing with a basic Sears- 
Haack body and a body indented symmetrically for a Mach number of 1.2 
is shown in figure 1. The other two bodies, one with a symmetrical 
M a 1.4 indentation and one with an asymmetrical M = 1.4 indentation, 
are shown together with the wing in figure 2. The symmetrical M = 1.4 
body has a cylindrical shape in the region of the wing-body juncture. 
Consideration of the effects of angle of attack led to the design of the 
asymmetrical body, which is indented on the upper surface but which 
actually has a bulge on the lower surface. Ordinates for all the bodies 
are presented in table I. 

For the remainder of this report, the four bodies will be desig- 
nated as follows: The basic Sears-Haack body will be called the basic 

body; the body indented symmetrically for a Mach number of 1.2 will be 
called the symmetrical M = 1-2 body; the body indented symmetrically 
for a Mach number of 1.4 will be called the symmetrical M = 1.4 body; 
and the body indented asymmetrically for a Mach number of 1.4 will be 
called the asymmetrical M = 1.4 body. 

For the transition- fixed tests, 0.10- inch transition strips were 
located at 10 percent of the chord on both upper and lower surfaces of 
the wing and at 10 percent of body length. The strips were obtained 
by spraying the surfaces with a commercial liquid plastic and blowing 
on carborundum grains (approximately 0.012 inch in diameter) at an 
estimated density of 40 grains per inch. 
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DATA AVAILABLE FOR ANALYSIS 


Force data have been published in reference 6 for the basic and 
1.2 symmetrically indented configurations with both fixed and nat- 
ural transition at Mach numbers from 0.80 to 1.45. Pressure distribu- 
tions have been published in reference 7 for these same two configura- 
tions with transition fixed. & 


Additional pressure distributions for the easic and M = 1.2 sym- 
metrical indentations with transition natural, md for the symmetrical 
and asymmetrical M = 1.4 indentations with transition fixed were 
obtained in recent tests at a Mach number of 1 . 45 . These tests were made 
in the Langley 8- foot transonic pressure tunnel where the supersonic Mach 
number of 1.45 is obtained by enclosing the longitudinal slots with spe- 
cially designed channels that convert the slotted test section to a 
supersonic nozzle. Details of the resulting nozzle shape and of the 
test-section Mach number distributions have been published in reference 8. 

, r these tests , the angles of attack were -2°, 0°, 2°. and 4 °. Pressure 
distributions were obtained at six wing stations. 12. 25 40 60 80 

and 95 percent of the wing semispan, and at five longitudinal rows of 
body orifices, which were spaced at intervals of 45° and designated as 
rows A, B C, D and E. The rows of body orifices were lettered starting 
ith row A on the center line along the upper surface and progressing in 
a clockwise direction while facing downstream. The orifice locations 
are given in reference 7. The results of these tests are presented 
herein in the form of pressure coefficients in figures 5 and 4 


METHOD OF ANALYSIS 


The purpose of this section is to describe the methods by which 
the wave drag due to thickness, the wave drag dua to lift, the vortex 
drag, the integrated pressure drag, and the skin -friction drag were 
either calculated, op derived . 


Wave drag .- Area distributions for the comb 
body and the symmetrical M = 1.2 body have bee: 
ence 6 for Mach numbers of 1 . 0 , 1 . 2 , and 1 . 4 . T 
obtained in the normal manner (ref. 2) are ident 
angles 0 and l80° - 9. Furthermore, they are 
butions that would be obtained if the wing were ; 
The wave drag due to thickness has been calculate 
tributions by means of the method of reference 9 


Lnations with the basic 
1 presented in refer- 
lese area distributions 
cal for cuts taken at 
the same area distri- 
symmetrical (uncambered) 
; d from these area dis- 


Calculation of the wave drag for a lifting wing-body combination 
is more complex than for a nonlifting combination. In reference 5, the 
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wave drag is shown to be a function of S"(x) - (f3Z ! /2q), where S M (x) 
is the second derivative of the cross-sectional area distribution and 
Pi'/2q is a function of the pressure distribution. 

The calculation of the wave drag then requires the determination 
of both the cross-sectional area distribution and the distribution 
along the streamwise axis of the quantity Z . For the 0 = 0p cut, 

this quantity Z is the component, normal to the free stream and paral- 
lel to the plane 0 = 0^, of the net resultant force on a section formed 

by the intersection of a Mach plane for the 0q cut with the configura- 
tion surface (ref. 5) * 

In the present analysis, several simplifying assumptions have been 
made in determining the distributions of J3Z/2q. These may be listed as 
follows : 

(1) The wind axes specified in the theory were replaced by body 
axes in the calculations for an angle of attack of 4°. 

(2) The side force determined by the integration of the pressures 
in the y-direction was neglected in obtaining the net resultant force. 
(For one case that was checked, calculations of the magnitude of this 
side force for 0 = 90 °, at which point it would be a maximum, indicated 
that the net effect on the drag was negligible.) 

(3) In determining the streamwise distribution of (3Z/2q, only the 
region between the wing apex and the intersection of the cutting plane 
from the wing-tip trailing edge with the body center line was considered 
(Calculations for a = 4° using a distribution of this type and also a 
distribution that started at the body nose showed that omitting the lift 
on the forebody had only a small effect on the wave drag.) 

Utilizing the second of these assumptions and the fact that the 
wing-body combination is symmetrical about the x,z-plane, the quantity 
Z 1 can be shown to have the following property: 


that is, Z 1 (x, 0) 
me try, the term 


X, 0 - 


7 

2 ) 

. rt 

e. 

2 ” 


S ,r (x,0) is seen to be an even function of - 0. If 


the terms S n (x,0) and Z*(x, 0) have the abovementioned properties, 
reference 5 shows that the wave drag due to thickness and the wave drag 
due to lift can be analyzed separately. All the wing-body combinations 
discussed in this paper will be analyzed in this manner. The error 
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in calculating the wave drag due to lift (introduced by the above men- 
tioned assumptions) has been determined and was found to be less than 
2 percent. 

The (3l/2q distributions , as determined from pressure distributions 
at an angle of attack of 4°, are shown in figure 5 for the wing-basic 
body combination at several Mach numbers and values of 0. Distributions 
of (3Z/2q for the asymmetrical M = 1.4 body and wing combination at a 
Mach number of 1.43 and for the symmetrical M = 1.2 body and wing com- 
bination at Mach numbers of 1.43 and 1.125 are compared in figure 6 with 
the distributions for the basic body. The incremental wave drag due to 
lift was calculated for each distribution by the method that is given in 
reference 9 for calculating the wave drag due to thickness from distri- 
butions of S'(x). The distribution of the ir.cremental wave drag due 
to lift with 0 was then integrated to obtair the total wave drag due 
to lift. 

Recently, it has been suggested (ref. 4) that for asymmetrical wing 
locations and for cambered wings the wing should be considered as a reflec- 
tion plane and the areas above and below the wing- chord plane should be 
treated separately. This procedure is based on the consideration of a 
wing with an attached leading-edge shock wave, so that there is no com- 
munication between the upper and lower surface flow fields. 

Let the area slope distribution for the 0 = 0]_ cut, by the method 
of reference 4, be given by 

S'(x,0i) t + S'(x,9 1 ) a 

where sl ( x jGi)t is 3. eriv ed from the area distribution of the complete 
model in the normal manner (ref. 2). The tern SI ( x j 0 i) a is the con- 
tribution of the configuration asymmetry about the wing-chord plane and 
is derived from twice the difference in the aiea above and below the 
wing- chord plane. 

Area distributions due to configuration e symmetry for a Mach num- 
ber of 1.4 for the four wing-body combinations are presented in figure 7- 
The wave drag due to the configuration asymmetry has been calculated on 
the basis of these area distributions and the results will be presented 
in this paper. 

It is now suggested that the method of reference 4 might be extended 
to the case of a wing-body combination at an angle of attack and that an 
equivalent area distribution due to angle of attack be found. This pro- 
cedure will also be based on the consideratior of a wing with an attached 
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leading-edge shock wave, so that there is no communication between the 
upper and lower surface flow fields . 


The wing is assumed to be a flat plate (in the wing-chord plane) 
with an angle-of-attack distribution over the semispan. This distribu- 
tion of angle of attack may be any combination of angle of attack of the 
body, wing incidence, wing twist, and effective twist due to body induced 
upwash. The cross-sectional area due to angle of attack for any value 
of x on the body axis is twice the projection normal to the free stream 
of the area on the cutting plane between the wing- chord plane and a plane 
parallel to the free- stream direction passing through the wing leading 
edge (see fig. 8). This area may be expressed mathematically as the 
integral of the height h on the projection over that portion of the 
semispan in which the cutting plane intersects the wing-chord plane or 


S(x) a 



h dy 


From figure 8, it can be shown that h = -£ tan a, and for small values 
of a, h = -£a. Nov 

£ = (cot t - tan A)y + x 


Thus, 


S(x) a 



tan A)y + xja dy 


Differentiating with respect to x yields 

>y 2 


i*(x)a - -2 / 

J y± 


a dy 


Therefore, the slope of the cross-sectional area distribution due to 
angle of attack of the wing is equal to twice the area under the curve 
of spanwise distribution of angle of attack between the limits y^ and 

Y 2 - These limits are the points of interception of the Mach cutting 

plane with the boundaries of the exposed wing- chord plane, the inboard 
point being the lower limit. 


The wave drag due to angle of attack has been calculated from the 
area slope of several of the wing-body combinations used in this investi- 
gation and the results have been presented in this paper. 
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Vortex drag .- The wave drag is derived in reference 5 by considering 
the transfer of momentum across a cylindrical surface of large radius that 
surrounds the lifting wing-body combination and extends far enough to the 
rear so that all the wave drag appears on the cylinder. The remaining 
drag can be found by considering the transfer of momentum across the end 
of this cylinder , a disk at a large distance downstream. At subsonic 
speeds, this part of the drag is normally referred to as induced drag, 
but this terminology seems inappropriate at supersonic speeds. The term 
vortex drag will be used throughout the entire speed range in this paper. 

If 1 the spanwise distribution of circulation I" 1 is known, the vortex 
drag can be evaluated by the following equation: 


c D,v 


1 

2 *U 2 S w 



pb/2 

J_ b/2 r ^ r '( y l> log ek - 


y-J dy dy L 


This equation has the same form as the wave-drag equation (ref. 5) and 
consequently , the methods of reference 9 were used in its solution. 


The spanwise distribution of circulation was found from the relation 


r = 


c av U 



where c n was determined from the pressure distribution over the 

u av 

surface of the wing-body combination (ref. 7). 

Integrated pressure drag .- The integrated pressure drag coefficients 
were calculated from the pressure distributions by use of the following 
equations: 



where c^ - c x cos a s + c n sin a g . The wing twisted appreciably under 

load because of aeroelastic effects (ref. 7). The section angle-of- 
attack distribution was calculated for each Mach number and angle of 
attack for the wing — basic-body combination using the section force 
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coefficients and wing-twist influence coefficients tabulated ^refer- 
ence 7 The section drag coefficients were then obtained by adding the 
comp onent s^ 6 in the drag direction, of the section chord-force and normal- 

force coefficients. 

Skin-friction drag- The skin- friction drag throughout the Mach num- 
— - - r tn \-— from the overall Cp 


ber range was obtained by subtracting ( c D,p)v/B 
obtained in the force tests of reference 6. 


RESULTS AND DISCUSSION 


Pressures Induced by Body Indentation 

The difference between the pressure on a wing— basic-body combination 
and thf pressure on a wing-indented-body combination can be considered 
as a pressure induced by changing the body shape. Contours of consta 
induced pressure coefficients dC p axe shown in figure 9 for the M . 1.2 
indented-body combination at 0° angle of attack and several Mach numbers . 
The contours were obtained from pressure data in reference 7- 

P ressures induced on body - The differential body pressures in fig- 
ure 9 show as would be expected, an expansion over the forward par 
the indentation, followed by a compression. The strength of the expan- 
sion on the body increases with increase in Mach number, becoming a 
maximum at M = 1.03, and then decreases with further increases in Mach 
number. The compression, however, is strongest at subsonic spee s an 
decreases with increase in Mach number. In both cases, the strength of 
the induced pressure field appears to be a maximum when the local Mac 
number in the vicinity of the body is approximately 1.0. 

In reference 11, the induced pressure field is calculated by lin- 
earized theory for indented and bumped bodies for which the source 
distribution's given. When the body shape, rather than the source 
distribution is given, the calculation of the pressure field is more 
aimcSS in thf limiting case on the surface of the body, however, 
the change in pressure coefficient can be obtained from e wo- 
dimensional relation (ref. ll) 


- f 


where Is the change in body slope. However, the pressure change 

calculated by this two-dimensional relation is displaced downstream. 
The measures induced over the upper surface of the body by the sym- 
StrSal T; 12 indentation are compared in figure 10 with induced 
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pressures given by the two-dimensional relation for Mach numbers of 1 12b 
and 1 .43. The calculated pressure distributions are not only displaced 
downstream, as expected, but the pressure coefficients also are of greater 
magnitude than the experimental values. Both theory and experiment show 
a decrease in induced pressure coefficient with Mach number. 

_. Pressures induced on wing.- The location of the induced pressure 
field relative to the wing plan form varies greatly with free- stream 
Mach number (fig. 9 ). At a Mach number of 0.9^, the induced pressure 
field extends more or less laterally across the wing. With increase in 
, C f r ^ rom the expansion field from the forward portion of 

e indentation increases in spanwise extent anu is swept back along the 
w ng until, at a Mach number of 1.43, it covers the major portion of the 
W J nfi - L ^ kewi se, the compression field, which i £ well forward on the wing 
at a Mach number of O.94, is gradually swept back with increase in Mach * 
number until it affects only a small portion of the inboard trailing- 
edge region of the wing at a Mach number of 1.43. The strength of the 
pressure fields over the wing is attenuated with spanwise distance. At 

WMCh 13 nearest the indentation design Mach num- 
ber of 12 , the induced pressure fields are felt at the wing tips, but 
due to the attentuation, the major effects are limited to the inboard 
60 percent of the semispan. 


ZNCp = 0 


Line of zero pressure change- The line of zero pressure change 
: 0 is of particular interest because it represents the boundary 

between the expansion and compression fields. Examination of the actual 
^efficients on the wing— indented-body combination indicates 
that this line is generally swept back so that it is tangent to the Mach 
lines associated with the local velocities over the wing and body. 

0 viously, the body indentation will be most effective in reducing the 
wmg drag if this line is swept so as to coincide with the line of maxi- 
mum wing thickness. At a Mach number of 1.125, -he ££ =0 line on 

the upper wing surface (fig. 9 ) approaches this optimum** location (maxi- 
mum thickness occurs at the 40-percent chord for this wing) . However 
because of differences in the flow fields due to camber over the upper 


lower surfaces of !h.e wing^ fhe 
is located forward of the Z£ p = o 


Z£ p = 0 


line on the lower surface 


line on the upper surface. 

ang le of attack .- Figure 11 shows contours of 
induced by the symmetrical M = 1.2 indentation for a = 4° angle of 
attack and Mach numbers of 1.125 and 1.43. Only small differences are 
noted in the strength of the pressure fields at C° and 4° angle of 
attack. The main effect is the further displacement of the ££ = 0 

lines, and consequently the pressure fields, on the upper and lower sur- 

faces of the wing. This displacement results in a small redistribution 
of the lift on the wing and body. 
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Effect of asymmetrical indentation .- The asymmetrical M = 1.4 
indentation was designed specifically to produce a redistribution of 
lift. The pressure fields induced on the upper and lower surfaces of 
the wing by this asymmetrical indentation are shown in figure 12 for 
4° angle of attack and a Mach number of 1.43 • The induced pressures 
shown are those which resulted from changing the symmetrical M = 1.4- 
indentation to an asymmetrical indentation. It should be mentioned 
that the total cross-section area distribution for the two configura- 
tions are different (see ref. 6). However , it is felt that the differ- 
ences are not of sufficient magnitude to obscure the effects of asym- 
metrical indentation. The asymmetrical indentation produced an expansion 
on the upper surface and a compression on the lower surface over the 
forward portion of the body and a compression on the upper surface and 
an expansion on the lower surface over the rearward portion of the body. 
The compression field on the lower surface of the asymmetrical indenta- 
tion is relatively weak and is actually ineffective in producing lift. 

The result of the asymmetry on the induced pressure fields is to increase 
the lift over the forward regions of the wing and to decrease it over the 
rearward regions. „ This redistribution of lift will be shown subsequently 
to have a favorable effect on the wave drag due to lift. 


Drag Analysis for a = 0° 

Spanwise distribution of drag .- The spanwise variations of section 
drag loading coefficient for the wing in the presence of the basic and 
symmetrical M = 1.2 bodies for several Mach numbers and for . the sym- 
metrical and asymmetrical M = 1.4 bodies for a Mach number of 1.43 are 
presented in figure 13 . The symmetrical M = 1.2 indentation reduced 
the drag loading near the wing-root juncture at all Mach numbers. At 
Mach numbers of 0.94 and 0.98 body indentation increased the drag of 
the outboard sections. This adverse effect was due to the forward 
location of the compression field at the low Mach numbers. With increase 
in Mach number the reduction in drag loading increased in spanwise extent 
as the induced pressure fields associated with the change in body shape 
were swept back across the wing. However, the major drag reduction gen- 
erally occurred over the inboard 60 percent of the semispan. The sym- 
metrical M = 1.4 body, which was cylindrical in the region of the wing, 
reduced the drag of the wing but was not so effective as the M = 1.2 
indentation. A further, but small, reduction in drag loading was 
obtained when the body was indented asymmetrically. The average drag 
loading over the bodies at subsonic speeds was lower than the drag 
loading of the wing near the wing-body juncture; but, it became con- 
siderably higher at supersonic speeds. 

All the data which have been discussed up to this point are for 
transition fixed on the wing and body. The pressure distribution over 
the wing and body for a Mach number of 1.43 with transition natural and 
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with transition fixed are compared in figure 3 for the basic and sym- 
metrical M = 1.2 bodies. Fixing transition on the wing generally 
increased the leading-edge pressure coefficient at 0 ° angle of attack. 
Figure 12 shows that this increase in leading-edge pressure caused a 
drag increase over the inboard regions of the wing. 

Integrated pressure drag for wing and body .- The pressure drag for 
the wing has been obtained by integrating the section drag loading 
between the wing-body juncture and the wing tip, and the drag for the 
body by performing the integration indicated in the section entitled 
"Method of Analysis." The variation of the wi:ig pressure drag, body 
pressure drag, wing-body pressure drag, and skin- friction drag with 
Mach number are presented in figure 14 for the basic and symmetrical 
M = 1.2 bodies with fixed transition. Body indentation had no effect 
on the drag of the wing for Mach numbers below 0.94. Above this Mach 
number, indentation was effective in reducing the drag at all Mach num- 
bers. The drag for the M = 1.2 body was higher at a Mach number 
of 0.80 than that for the basic body. The drag-rise Mach number for 

the indented body was delayed and the drag rise reduced at all Mach 

numbers except in the region of 1.43. For the wing-body combination 

body indentation increased the drag at subsoni : speeds and reduced the 
drag rise at all Mach numbers tested. The values of the wing pressure 
drag, the body pressure drag, and the wing-bod;/ drag for a Mach number 
of 1.43 are listed in the following table for all four configurations 
with transition fixed and for the basic and symmetrical M = 1.2 wing- 
body combinations with transition natural: 



Transition fixed 

Transition natural 


Basic 

body 

Symmetrical 
M = 1.2 
body 

Symmetrical 
M = 1. 14- 
body 

Asymmetrical 
M = 1.4 
body 

Basic 

body 

Symmetrical 
M = 1.2 
body 

c D,p,W 

0.0077 

0.0056 

0.0070 

O.OO 65 

0.0071 

0.0051 

C D,p,B 

• 0035 

.0036 

.0031 

.00^7 

.0028 

.0031 

C D,p,WB 

.0112 

.0092 

.0101 

. 01 L 2 

.0099 

.0082 

C D,f 

.0066 

.0078 


— 

.0071 

.0063 


For the asymmetrical M = 1.4 body, the irag of the wing, as dis- 
cussed previously, is reduced. However, it cai be seen from the pre- 
ceding table, that the drag of the body is higier than that of any of 
the other bodies. As a result, the asymmetrical M = 1.4 wing-body 
combination has a drag which is higher than that of the symmetrical 
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M - 1 4 wing-body combination and is equal to that of the basic wing- 
body combination. The effect of transition on the distribution of wing 
section drag at a Mach number of 1.45 has been discussed previously. 

With fixed transition, the drag for the wing increased for both the 
basic and symmetrical M * 1.2 wing-body combinations. The drag for 
the bodies increased also (see preceding table). The net result is 
increase in drag of approximately 15 percent. 

Skin- friction drag.- The variation of skin- friction drag with Mach 
numbe r (fig. l4) for ~bfth the basic and symmetrical M - 1.2 wing-body 
combinations with transition fixed was obtained from Pressure data 

obtained in this investigation and the force test results of reference 6. 
The skin- friction drag for the basic wing-body combination decreased 
sharply with increase in Mach number at transonic speeds and increased 
slightly at supersonic speeds. The drag for the wing— indented-body 
combination was lower at low subsonic speeds than that for the wing- 
basic-body combination; but, with increase in speed the skin-friction 
drag for the indented wing-body combination was higher. These 
SboS shape on the friction drag, although generally small must be 
assumed to be due to the influence of the induced pressure fields on 
the boundary layer. At a Mach number of 1.45, fixing transition had 
only a small effect on the wing— basic-body combination but produced a 
SSe incase in skin-friction drag for the wing- indented-body com- 
bination (see table given in preceding section). It is evident on the 
basis of these results that in obtaining the pressure drag rise fro 
force data the assumption of a constant skin-friction drag at transonic 
speeds evk for transition fixed, could in some cases lead to erroneous 
results. For the two wing-body combinations being considere , ® 

rise obtained in this manner from the force data indicates no effe 
IX indentation at a Mach number of 1.45 (ref. 6); actually the drag 
rl7e obtained from pressure data indicates that body indentation pro- 
diced a considerable reduction in wave drag at this Mach number For 
this reason, the calculated drag values that are presented in this 
paper are, wherever possible, compared with drag values obtained from 

pressure data. 

Comoarison of calculate d drags with integrated pressure dra g- A 
comparison of the calculated drags with the integrated pressure drags 
fo/the basic and symmetrical M - 1.2 body and wing combinations at 
an angle of attack of 0° is presented in figure 15- The calcula 
drags consist of the vortex drag, the wave drag due to thickness, and 
the wave drag due to lift. Values of these various drag components were 
obtained in the manner described in the section entitled Method of 
Analysis ." 

For both wing-body combinations at subsonic speeds (M - °-9°) th ® 
calculated vortex drag underestimates the integrated pressure ag. 
transonic speeds (0.94 ^ M < 1.125) the sum of the calculated drags does 
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combinSion ely TMf iCt . integrated pressure drag for either wing-body 
. , r; ^ * Thls result is not suprising since both the wave drag due 

to thickness and the wave drag due to lift were calculated by means of 

soni“Jeloclty er C WhlCh d0eS n0t app1 ' at speeds cl ° se 

At the supersonic Mach number of 1 . 43 , the sum of the calculated 

thp 6 ? +° r t + e basic bod ^ and win e combination is only 0.0004 lower than 
he integrated pressure drag. For the same Mach number the sum of the 

if 0 U 0015 lowef thf S® ay r etrical body and wing combination 

IS 0.0015 lower than the integrated pressure drag. A large portion of 

this error might well be attributed to the inadequate prediction of the 

dT ^ g ’ ^ iCh 13 evident at the subsonic Mach numbers and probably 
extends throughout the Mach number range. ^ ^ 

In the section entitled "Method of Analysis' it was stated that 
the supersonic wave drag is a function of the axial distribution T 
tv \<l /2q} ' F a s y mmetric al nonlifting wing-body combination 

calculated and included in the comparisons of figure 15. For both the 

ffonf? S f72f 1 afoO M = f 2 " 0dy dofflbinations^the^dis tri- 
presen ted in figure 16 f“ 8 9 f 0° anfiff S number ° f ^ 1= 

The effect of symmetrically indenting the basic bodv and uinn 
bination for a design Mach number of 1.2 is shown in figure 15(c ) & The 

sSefiff ' M 5 Tl “"f f pressure ^ag between the bash and 
Cai M “ 1,2 Icdy and wing combinations. The solid lina rmw 
sents the difference between the total calculated drags of the two com 

witfrf D . 0 . 4 3 of tfzz . 

accuracy throughout the ifcfnufe f ranged omf .ftffS “f efaf thf^ 
d“f 4 fcfacy f 

Approximate method for determining wave rj-rncr ^ . 

S BUi- bhe suggestion is made in fefn e T St t^Tf/ecTj^ 
Wing camber can be taken into account in the supersonic area rf e f con 
sidering the wing-chord plane as a reflection plan 2 and treating seL 

Jr : n f s s f ssss f ^^fffhffieLf 
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the wing-body combinations of this paper since, even at M = 1.43, the 
wing leading edge is swept behind the Mach line. However, in reference 12 
it is shown that wings with subsonic leading edges can exert a powerful 
effect as a divider plate. Therefore, it appears instructive to deter- 
mine the value of this method for predicting the increment in wave drag 
due to asy mme try about the wing-chord plane. The question arises as to 
whether this method will give the same result in practice as would 
accounting for the camber by including the PZ/2q term. 

A comparison of the slope of the area distribution due to camber 
(as determined by the method of ref. 4) and -PZ/2q is presented in 
figure 17. Also shown is a modification of the method of reference 4. 

The modification is obtained by multiplying the slope of the area dis- 
tribution due to camber by cos 0. This step makes the analogy between 
the method of reference 4 and the theory of reference 5 more complete 
since the lift distributions pz/2q are multiplied by cos 0 in order 
to obtain the component of force in the 0-direction. In determining 
S ' ( x) camber an a ll° warice h as been made for a small amount of aero- 
elastic^wing twist by an approximate method that will be discussed 
subsequently . 


Comparison of the curves (fig. 17) indicates poor agreement between 
S'( x ) camber and P^/ 2c 1* The agreement is not helped by multiplying 

S*( x ) camber by cos 0 - If the flow fleld over the wing Were two_ 
dimensional, that is, if the pressure coefficient at any point were pro- 
portional to the slope of the surface at that point, the distribution of 
pz/2q would be correctly predicted by S ' (x) cara ^ er cos 0. The flow 

field is not two-dimensional, however, since the wing is swept behind 
the leading-edge Mach line and the wing lies within the region of influ- 
ence of the body . 

The wave drag due to camber has been calculated for the wing — basic- 
body combination at M = 1.43 and a = 0° from the lift distributions 
(pz/2q) and from the modification of the method of reference 4 
(S ' (x) camber cos e )- The resultant drag-coefficient increments are 

0.00108 and 0.00025, respectively. 

In the case of the asymmetrical M = 1.4 body, the drag analysis 
requires an estimate of the contribution of the body asymmetry to the 
wave drag at a = 0°. If the lift distribution is known, the incremental 
drag can be calculated from the distributions of pz/2q. The suggestion 
is made in reference 4 that the effect of body asymmetry can be accounted 
for in a manner similar to that used for wing camber. 

A comparison of the slope of the area distribution due to wing cam- 
ber and body asymmetry (as determined by the method of ref. 4) and the 
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distribution of (3Z/2q. is presented in figure 18. Also shown is a 
modification of the method of reference 4- obtained by multiplying the 
slope of the area distribution due to wing camber and body asymmetry by 
cos 0 . The reason for this modification was discussed previously. 

Figure 18 indicates that neither the methci of reference 4 nor the 
modified method can predict, even approximately, the distribution of 
3 Z/ 2 q for a wing-body configuration on which the area is distributed 
asymmetrically about the wing-chord plane. The values of wave drag due 
to lift calculated from the distribution of 3 z/ 2 q and the modified 
method are O.OOO 58 and 0.00047, respectively. 


Drag Analysis for a = 4 3 

Spanwise distribution of drag .- The chordwLse pressure distributions 
presented in reference 7 have been integrated aLong the wing-section 
thickness to obtain the section chord force. Tils chord force, and the 
section normal forces listed in reference 7 were resolved about the wind 
axes to obtain a section drag coefficient. Because of the flexibility 
of this wing, appreciable wing twist due to aeroelastic effects occurred 
at low angles of attack. The spanwise distribution of wing twist was 
calculated using the influence coefficients and aerodynamic loads pre- 
sented in reference 7 . These calculated twist angles were then used in 
resolving the forces about the wind axes. 

The spanwise variation of the wing-section drag loading coefficient, 
c d is sh o w n in figure 19 for Mach numbers of 1.03, 1.125, and 1.43, 

for the wing in the presence of a basic and a symmetrical M = 1.2 
indented body. Similar distributions are shown for the symmetrical and 
asymmetrical M = 1.4 indented bodies at a Mach number of 1.43. The 
symmetrical M = 1.2 indentation reduced the drag loading over the wing 
at Mach numbers of 1.03 and I.I 25 . At a Mach number of 1.43, however, 
the drag loading near the wing-body juncture is reduced; but, since the 
pressure field induced by the indentation is swept back sharply, the drag 
loadings for the outboard wing sections are increased. Both the sym- 
metrical and asymmetrical M = 1.4 indentations redistributed the drag 
loading over the wing, and caused relatively lai ge increases over the 
midspan. The drag for the asymmetrical indental ion was lower over the 
inboard 75 percent of the wing semi span than it was for the symmetrical 
indentation. 

Integrated pressure drag for wing and body .- At subsonic speeds, the 
body drag was increased when the basic body was indented symmetrically 
for a Mach number of 1.2 (fig. 20); however, sulstantial reductions 
occurred at supersonic speeds. The reduction in drag on the body due 
to the indentation generally was greater than that for 0 ° angle of attack 



19 


and can be attributed to a favorable interaction of the pressures induced 
on the body by the wing and the body slopes. For the asymmetrical inden- 
tation, this interaction was unfavorable and the drag of the body 
increased as shown in the following table for M = l.Vj: 



Basic 

body 

Symmetrical 
M = 1.2 
body 

Symmetrical 
M = 1.4 
body 

Asymmetrical 
M = 1.4 
body 

C D,p,W 

0.0180 

0.0177 

0.0176 

0.0164 

C D,P,B 

.0055 

.0050 

.0056 

.0066 

C D,p,WB 

.0235 

.0227 

.0232 

.0230 


The drag of the wing was reduced at all Mach numbers when the wing 
was in the presence of the symmetrical M = 1.2 indented body. However, 
because the £C p = 0 line is swept back farther along the wing at angle 

of attack, the reduction in wing drag at 4° angle of attack was con- 
siderably less than it was for 0° angle of attack. For the asymmetrical 
M = l.lt- indentation, the C p = 0 line is located farther forward on 
the wing and the preceding table shows that the largest drag reduction 
for the wing at a Mach number of 1.45 occurred for the wing in the pres- 
ence of the asymmetrical M = 1.4 indented body. 

Because of the reduction in effectiveness of body indentation on 
the wing drag, the total drag reduction due to indentation was generally 
smaller at 4° angle of attack than it was at 0° angle of attack (figs. 14 
nnrt 20). The drags for the basic and symmetrical M = 1.2 combinations 
as obtained from the pressure data also are compared in figure 20 with 
the force test results of reference 6. Because of differences in rigidity 
of the two wings, the incremental drag between the two sets of data 
included not only skin- friction drag but also drag differences due to 
aeroelastic wing twist . 

Vortex drag.— The spanwise load distributions at several Mach num- 
bers for the wing-body combinations with the basic and symmetrical 
M = 1.2 bodies are presented in figure 21. Also shown are the span- 
wise load distributions for the symmetrical and asymmetrical M = 1.4 
indented body and wing combinations at a Mach number of 1.43- The vor- 
tex drag has been calculated from these distributions . The results are 
presented in figure 22 in terms of the vortex drag due to lift parameter 

C p v j . The variation with Mach number is much the same for both the 

basic and symmetrical M = 1.2 body and wing combinations. The value 
of the vortex drag due to lift parameter varied from 42 to 55 percent 
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above the value of 1/itA, the value for an elliptic span-wise load dis- 
tribution. The value of the vortex drag due vo lift parameter for the 
symmetrical M = 1.4 and asymmetrical M = 1.4 indented- body and wing 
combinations are both about 51 percent above '.he minimum value of l/nA. 

Wave drag due to lift .- The method used in calculating the wave drag 
due to lift has been dsicussed in the section entitled "Method of 
Analysis." It is believed, however, that the distributions of pi/2q 
presented in this paper are perhaps the first experimental distributions 
obtained for a sweptback wing and warrant considerable discussion. 

The required distribution of pz/2q for the minimum drag due to 
lift derived from the wave-drag equation of reference 5 is elliptical 
for all values of 9. A similar conclusion we s arrived at by Jones in 
deriving the minimum drag due to lift for the special case of an 
e lliptl ca l - pl an_ f'° rm wing (ref. 15). The experimental distributions 
obtained for the 45° sweptback wing — basic-body combination depart 
greatly from elliptical shapes and vary considerably with 9 and with 
Mach number (fig. 6). For a Mach number of 1.45, the 9=0° lift 
distribution is smooth but builds up rapidly to a sharp peak. With 
increasing 9, the distributions are elongated along the body axis; how- 
ever, the peaks of the distributions tend to tecome more abrupt. The 
magnitude of the peaks generally decreases with increasing 9 because 
the Z/ q term is a function of the cosine of the angle 9. With decrease 
in Mach n umber, the magnitude of the distributions decreases because of 

the P = |/m 2 - 1 term. The actual shapes of the distributions also 

t)6C£ius6 of the effect of Ms-ch number on the cutting plane angle . 

For 9 = 0°, the only effect is in the p tern; and, the shapes of the 
distributions for a constant lift are essentially the same. 

The incremental wave drag calculated from the 9-distributions for 
several body shapes are plotted against 9 in figure 25 for Mach numbers 
of 1.45 and 1.125. The symmetrical M = 1.2 indentation, as noted pre- 
viously, caused a small redistribution of the Lift on the wing at angles 
of attack. Thus, although the wave drag due ti lift is independent of 
the thickness drag, the redistribution of lift caused by indenting the 
body to obtain a reduction in thickness wave d~ag can influence the wave 
drag due to lift. Indentation, at a Mach number of 1 . 125 , caused a 
small reduction in the peaks and a small incre ise over the forward por— 
lion of the P Z/2q distributions for the wing —basic— body combination 
(fig- 5) • These changes in the pz/2q distributions were favorable 
and, as indicated in figure 25, reduced the wave drag due to lift. With 
increase in Mach number to 1.45, the induced pressure fields on the upper 
surface of the wing are swept back sharply (see fig. 10) and the addi- 
tional lift which results tends to aggravate tie peaks of the pz/2q 
distributions and increase the drag. The asymmetrical M = 1.4 inden- 
tation was designed specifically to alter the lift distribution of the 



21 


wing. The peaks of the PZ/2q distributions for both 0 cuts were 
reduced by the asymmetrical indentation and a substantial drag reduc- 
tion occurred. 

The parameter for wave drag due to lift v /c for the wing — 

basic-body combination is compared in the following table with the param- 
eter calculated for the symmetrical M = 1.2 and the asymmetrical 
M = 1.4 indented body and wing combinations. 


1 

Mach 

number 

Basic 

body 

Symmetrical 
M = 1.2 
body 

Asymmetrical 
M = 1.4 
body 

Lower 

bound 

1.03 

0.0089 

0.0087 


0.00473 ■ 

1.125 

.0413 

.0381 


.0152 

1 A3 

.1590 

.1454 

0.1134 

.0588 


Also shown in the table is the parameter for wave drag due to lift cal- 
culated for a wing-body combination which would have distributions of 
the same length as the experimental wing body but which had elliptical 
shapes. The drag for this hypothetical wing-body combination may be 
considered as representing a lower bound for the wave drag due to lift. 

It should be kept In mind that the minimum wave drag that can be achieved 
for a particular plan form wing and body combination may be higher than 
the lower bound. 

Comparison of calculated drags with integrated pressure drag .- A 
comparison of the sums of the calculated drags with the Integrated pres- 
sure drags for the basic and symmetrical M = 1.2 body and wing com- 
binations at an angle of attack of 4° is presented in figure 24. The 
calculated drags are the vortex drag, the wave drag due to thickness, 
and the wave drag due to lift. Each of these drag components has been 
discussed previously. 

For both wing-body combinations at subsonic speeds (M = 0.94) the 
vortex drag computed from the spanwise load distributions underestimates 
the integrated pressure drag by 6 to 13 percent. At transonic speeds 
(0.94 < M ^ 1.125) the sum of the calculated drags does not adequately 
predict the integrated pressure drag for either wing-body combination. 

As was the case for an angle of attack of 0°, agreement was not expected 
since both the wave drag due to thickness and the wave drag due to lift 
were calculated by means of linearized supersonic theory which does not 
apply at speeds close to sonic velocity. 
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At the supersonic Mach number of 1.43, the sum of the calculated 
drags for the basic body and wing combination is approximately 8 per- 
cent higher than that for the integrated pressure drag. For this same 
Mach number , the sum of the calculated drags fcr the symmetrical 
M = 1.2 body and wing combination is less than 1 percent below the 
integrated pressure drag. The sum of the calculated drags for the asym- 
metrical M = 1.4 body and wing combination is 0.0214 compared with the 
integrated pressure drag which is 0.0230. The sum of the calculated 
drags underestimated the integrated pressure drag by approximately 
7 percent. 

Approximate method for determining distribution of 3Z/2q .- In 
reference 4 the possibility of obtaining wave drag due to wing camber 
and body asymmetry has been suggested. The method has been modified in 
this paper in order to obtain a closer analogy between this method and 
the theory of reference 5* It is now suggested that the modified method 
might be extended to the case of a wing-body combination at an angle of 
attack and that an equivalent area distribution due to angle of attack 
can be found. The angle of attack may be any combination of angle of 
attack of the body, wing incidence, wing twist, and effective twist due 
to body induced upwash. The derivation of the equivalent area distribu- 
tion due to angle of attack is presented in the section entitled "Method 
of Analysis." 

The slope distributions of the areas due to angle of attack are 
plotted with the distributions of (3l/2q for the wing — basic-body com- 
bination at several Mach numbers in figure 24. The slopes of the areas 
due to angle of attack have been multiplied by 3 cos 0 in order to 
obtain a better comparison with variation in Mach number. These distri- 
butions include the effect of aeroelastic wing twist (ref. 7), effective 

/ 

wing twist due to body induced upwash (Oq = Org — as predicted in 

\ V y 2 

ref. 14 y, and wing camber (in the manner described previously in this 
paper) . 

At a Mach number of 1.03, the distributior s of pl/2q are predicted, 
within reasonable limits, by the negative slopes of the area distribu- 
tions due to angle of attack -3S ! (x)cos 0. At M = 1.125 the agreement 
between 3Z/2q and -3S f (x)cos 0 is quite gocd for 0=0°. For 
0 = 45°, the agreement is generally good except at the peaks (x ^ 12) 
of the distributions. For M = 1.43, the agreement for both 
0=0° and 45° is not too good, particularly toward the rear of the 
distributions. However, the magnitudes of the peaks and the general 
shape of the distributions of 3Z/2q are well represented by 
-PS'(x)cos 0. 
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Comparisons of pl/2q and -pS'(x)cos 0 for the symmetrical 
M = 1.2 and the asymmetrical M = 1.4 body and wing combinations at 
f/j = 1.45 are presented in figure 2 6 . Similar to the basic body and 
wing combination, the agreement is generally poor, particularly toward 
the rear of the distributions. Again, however, the general shapes of 
the Pl/2q distributions are well represented by -pS’(x)cos 9. 

Values of the incremental wave drag due to lift C D ^ Wj 0 have been 

computed from the distributions of pz/2q. and -pS'(x)cos 9 in fig- 
ures 25 and 26 and are presented in figure 27 as a function of 9 . The 
values of wave drag due to lift calculated from -pS'(x)cos 9 generally 
underestimate the values obtained from pi/ 2q. The differences were all 
within about 7 percent with three exceptions. Two exceptions occurred 
for the basic body and wing combination, 9 = 45° cuts, at 
M = 1.03 and 1.125. These differences were of the order of 0.0005 
and 0 . 0019 , respectively, which represents only 2 and 8 percent of the 
total calculated drag coefficients presented in figure 23- The third 
exception occurred for the symmetrical M = 1.2 body and wing combina- 
tion, 9=0° cut, at M = 1.43. This difference was 0.0041 and repre- 
sents l8 percent of the total calculated drag coefficient shown in 
figure 23- 

Values of the wave drag due to lift obtained from the curves of 
figure 26 are presented in the following table. Included in the table 
are the percent differences between Cp w calculated from -pS'(x)cos 9 

and Op w calculated from pz/2q (based on C^ w from pz/2q). These 
differences are also presented as percent of total calculated drag. 


Configuration 

M 

^D,w 

— 

Percent 

difference 

Percent 
of total 
calculated 
drag 

From PZ/2q 

From -pS'(x)cos 0 

Basic body 

1.03 

0.0011 

0.0008 

27.2 

1.3 

Basic body 

1.125 

.0042 

-0035 

16.7 

2.8 

Basic body 

1.43 

.0090 

.0083 

7.8 

2.8 

Symmetrical 

1.43 

.0093 

.0080 

14.0 

5.8 

M = 1.2 






body 






Asymmetrical 

1.43 

.0068 

.0066 

2.9 

.9 

M = 1.4 






body 
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The results in this table indicate that al'.hough the percent dif- 
ference between the wave drags due to lift calculated from -0S'(x)cos 9 
and 02/ 2q may be large, the errors introduced into the total calculated 
drag are reasonably small. For the five examples cited herein, these 
errors are all less than 6 percent. Thus, it appears that the method 
presented herein of obtaining the wave drag due to lift from the distri- 
bution of -0S'(x)cos 0 will afford an estimate of the actual wave drag 
due to lift without a knowledge of the pressure distribution on the 
configuration . 


CONCLUSIONS 


An investigation of effects of body shape on the drag of a 45° swept- 
back wing-body combination at Mach numbers from 0.90 to 1.45 resulted in 
the following conclusions : 

1. Both the expansion and compression fields induced by body indenta- 
tion were swept back as the stream Mach number increased from 0.94. The 
line of zero pressure change was generally tangent to the Mach lines 
associated with the local velocities over the wing and body. 

2. The strength of the induced pressure fields over the wing were 
attenuated with spanwise distance and the major effects were limited to 
the inboard 60 percent of the wing semi span. 

5- Asymmetrical body indentation tended to increase the lift on the 
forward portion of the wing and reduce the lift on the rearward portion. 
This redistribution of lift had a favorable effect on the wave drag due 
to lift. 

4 . Symmetrical body indentation reduced the drag loading near the 
wing-body juncture at all Mach numbers. The reduction in drag loading 
increased in spanwise extent as the Mach number increased and the line 
of zero induced pressure became more nearly alined with the line of 
maximum wing thickness. However, the major drag reduction generally 
occurred over the inboard 60 percent of the semispan. 

5- Calculations of the wave drag due to thickness, the wave drag 
due to lift, and the vortex drag of the basic ani symmetrical M = 1.2 
body and wing combinations at an angle of attack of 0° predicted the 
effects of indentation within 11 percent of the ving-basic body drag 
throughout the Mach number range from 1.0 to 1.45. However, the absolute 
values of the pressure drag were not predicted with the same degree of 
accuracy at Mach numbers of 1.125 and below. 
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6. Neither the method based upon the reflection of disturbances by 
the wing nor a modification of the method suggested in this paper suc- 
cessfully predicted the distribution of the oblique-section lift param- 
eter pi/ 2q for a wing-body configuration on which the area is distri- 
buted asymmetrically about the wing-chord plane. Consequently, these 
methods appear to be of little value in predicting the incremental wave 
drag due to lift resulting from wing camber and/or body asymmetry for a 
wing with subsonic leading edges . 

7. Calculations of the wave drag due to thickness, the wave drag 
due to lift, and the vortex drag for the basic, symmetrical M = 1.2, 
and asymmetrical M = 1.4 body and wing combinations predicted the 
total pressure drag to within 8 percent of the experimental value at 

M = 1.45. At Mach numbers of 1.125 and 1.03, the predictions were not so 
accurate. The incremental pressure drag due to a symmetrical M = 1.2 
indentation was predicted to within 8 percent of the total pressure drag 
of the wing — indented-body combination for all Mach numbers investigated 
above 1.0. 

8. A modification and extension of the method based upon the reflec- 
tion of disturbances by the wing to the case of a wing-body combination 
at an angle of attack of k° proved useful in predicting the distributions 
of pi/2q. The differences in wave drag due to lift calculated from the 
predicted and experimental distributions of pi/2q. were less than 6 per- 
cent of the total pressure drag. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., Aug. 1, 1958- 
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TABLE I.- Concluded 
BODY ORDINATES 


Body 
station, 
x, in* 

Radius, in* 

Basic 

body 

Symmetrical 

M » 1.2 
body 

_ 

Symmetrical 
M = l.k 
body 
(upper) 

Asym- 
metrical 
M = l.L 
body 
(upper) 

Asym- 
metrical 
M = l.h 
body 
(lower) 

18.5 

1.598 

l . lk 9 

1.586 

1.388 

1.689 

19.0 

1.602 

1.133 

1.586 

1.390 

1.69 k 

19.5 

1.606 

1.126 

1.586 

l.kOO 

1.697 

20.0 

1.606 

1.133 

1.586 

l . k 27 

1.701 

20.5 

1.60 k 

1.150 

1.586 

1.1*69 

1.702 

21.0 

1.602 

1.175 

1.586 

1.507 

1.700 

21.5 

1 1.600 

1.202 

1.586 

1.536 

1.692 

22.0 

1.59 k 

1.236 

1.596 

1.555 

1.675 

22.5 

1.587 

1.269 

1.591 

1.568 

1.6 k 9 

23.0 

1.578 

1.306 

1.576 

1.573 

1.622 

23.5 

1.570 

1.3 kl 

1.569 

1.568 

1.597 

2 k . 0 

1.560 

1.363 

1.55 U 

1.555 

1.57 k 

2 U .5 

1.5 k 7 

1.375 

1.539 

1.539 

1.5 k 9 

25.0 

1.532 

1.380 

1.521 

1.521 

1.521 

25.5 

1.517 

1.380 

i .5 oo 

1.500 

1.500 

26.0 

1.501 

1.376 

l . k 79 

l . k 79 

l . k 79 

26.5 

1.1*80 

1.370 

1.1*56 

1.1*56 

l . k 56 

27.0 

l . k 60 

1.362 

1.1*33 

l . k 33 

1.1*33 

27.5 ! 

1.1*38 

1.3 k 9 

l.klO 

1.1*10 

1.1*10 

26.0 

l.klk 

1.335 

1.388 

1.388 

1.388 

28.5 

1.387 

1.319 

1.365 

1.365 

1.365 

29.0 

1.360 

1.300 

1.3 k 0 

1.3 k 0 

1.3 k 0 

29.5 

1.330 

1.280 

1.315 

1.315 

1.315 

30.0 

1.300 

1.255 

1.290 

1.290 

1.290 

30.5 

1.269 

1.230 

1.259 

1.259 

1.259 

31.0 | 

1.231 

1.201 

1.229 

1.229 

1.229 

31.5 

1.193 

1.170 

1.195 

1.195 

1.195 

32.0 

1.158 

1.138 

1.156 

1.156 

1.156 

32.5 

1.116 

1.100 

1.116 

1.116 

1.116 

33.0 

1.076 

1.065 

1.073 

1.073 

1.073 

33.5 

1.030 

1.023 

1.030 

1.030 

1.030 

3 k . 0 

.98 k 

.980 

.985 

.985 

.985 

3 k . 5 

.936 

.933 

.93 k 

.93 k 

.93 k 

35.0 

.878 

.878 

.879 

.879 

.879 

35.5 

.825 

.825 

.823 

.823 

.823 

36.0 

.762 

.762 

.767 

.767 

.767 

36.15 

.750 

.750 

.702 

.702 

.902 






Figure 2.- Sketch of the wing with a symmetrical M = l.k and an asymmetrical 

indentation . 
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Figure 3.- Continued. 
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Figure 4.- Concluded. 
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Figure 7.- Distribution of cross-sectional area due to wing camber and 

body asymmetry at M = 1.4. 





(d) M = 1.125. 
Figure 9 .- Continued. 




(f) M a 1.43. 


Figure 9*- Concluded. 


Figure 10.- Comparison of the theoretical and experimental pressure 
change over the upper surface of the body due to a symmetrical 
M = 1.2 indentation. 
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(d) M = 1.43; lover wing surface. 


Figure 11.- Concluded. 
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(a) M = 1.43; upper wing surface 


(b) M = 1.43; lower wing surface. 

FI mire 12 - The pressure field induced over the wing and body by a 

^har^e' from a symmetrical M - 1.4 indentation to an asymmetrical 
indentation, a = b . 
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(b) 0 = 45° 


Figure 17 »- Comparison of the distribution of -3Z/2q and the slope 
of the cross-sectional area distribution due to wing camber, 
including the effects of aeroelastic wing twist, for the wing — 
basic-body combination, a = 0°; M = 1.43. 






Figure 18.- Comparison of the distribution of -pi/2q and the slope of 
the cross-sectional area distribution due to wing camber and body 
asymmetry, including the effects of aeroelastic wing twist, for the 
asymmetrical M = 1.4 body and wing combination, a = 0 ; M = 1.43. 
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Figure 20.- The variation with Mach number of the wing pressure drag, 
the body pressure drag, the wing-body pressure drag, and the force 
test results of reference 6 for the basic and symmetrical M = 1.2 
body and wing combinations . a = 4° . 
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Figure 22.- Variation with Mach number of the vortex drag due to lift parameter for the basi 
and symmetrical M = 1.2 body and wing combinations, a = 4°. 
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(b) Symmetrical M = 1.2 body. 

Figure 24.- Comparison of the experimental and calculated drags for the 
basic and symmetrical M = 1.2 body and ving combinations, a = 4°. 
Symbols indicate experimental values. 
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Fi glare 26.- Comparison of the distribution of pz/2q. and the slope of 
the area distribution due to angle of attack, wing camber, and body 
asymmetry for the symmetrical M = 1.2 and asymmetrical M = 1.1+ 
body and wing combinations . a = 1+ ; M = 1 .1+3 • 
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NASA - Langley Field, Va. 


to angle of attack, wing camber, and body asymmetry . 





